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ABSTRACT 

Doty, William, K., An Introduction to the Applicetion of Linear 
Programming to Industrial Engineering, May, 1954. There are 50 pages, 
5 figures, and 10 references in the bibliography. 

The purpose of this paper is to explain what linear programming 
is, how the linear programming problem is interpreted mathematically, 
how to solve a simple problem, what the assumptions, limitations, and 
possibilities of this technique are, and how to set up for solution a 
few examples of feasible problems. The emphasis is on Industrial 
Engineering applications and requires no knowledge of mathematics be- 
yond college algebra. 

The mathematical discussion illustrates the necessary elementary 
vector algebra including addition of vectors, unit vectors, and linearly 
independent sets. Linear transformations are illustrated by transfor- 
mations from a 3-dimensional to a 2-dimensional space. Convex sets 
and convex polyhedral cones are illustrated. A simple numerical example 
with three unknowns in two equations,with a linear functional to be 
maximized, is solved step by step and the geometrical significance of 
each step is explained. 

The ee necessary for the application of linear 
programming are discussed and certain feasible problem areas in 
Industrial Engineering are listed. It is shown why the technique can- 
not be used for scheduling, sequence of production selection, elimina- 
tion of expediters, etc. 3 

Three mathematical EE ace presented: (1) machine assignment 


to maximize total productivity, (2) production line loading for a single 








product to meet sales requirements over successive time periods and to 
minimize tho combined costs of regular and overtime production and 
storage, and (3) a numerical example to show the complexity to which 
even problems of small size are subject. This model deals with three 
products, two production operations, three time periods, and the max- 
imizetion of profit, 

The conclusion points out that the application of linear progran- 
ming to Industrial Engineering is, at present, limited, but predicts 
that, as study in both theory and practical aspects continues, more 


and more Industrial Engineering application will be made in the future. 





AN INTRODUCTION TO THE APPLICATION OF 


LINEAR PROGRAMMING TO INDUSTRIAL ENGINEERING 


INTRODUCTION 

Considerable interest has developed in the past few years in 
the mathematical technique known as linear programming. This interest, 
at first, centered itself among mathematicians and economists because 
the concepts were similar to or were advancements in technique under 
study in those two sciences. Interest in other fields was limitec 
until Dr. George B. Dantzig? evolved the "simplex method" of computa- 
tion. This was a turning point. The formerly theoretical and limited 
aspects of linear programming were now found to have direct application 


to practical problems facing the armed ón industry agriculture* 


and ,rensrportation^, Evidence of widening interest is illustrated by 
the appeerance of articles on linear programming in such practical 


6 
minded publications as Factory Management & Maintenance and Business 


en en ern a 


en nn nennen nn 


l. Dantzig, G.B., 'Maximization of a Linear Function of Variables 
Subject to Linear Inequalities," Activity Analysis of Production 
and Allocation, Cowles Commission Monograph 13; New York, John Wiley 
& Sons, 1951. pp 339-347. | | 

2. Wood, M.X., and Geisler, M.A., "Development of Dynamic Models for 
Program Planning," (ibid.), pp. 189-215. 

3, Charnes, A., Cooper, W.W., and Mellon, B., "Blencing Aviation 
Gasolines," Econometrics, vol. 20, r. 2. 

4. Hildreth, C., and Reiter, S., "On the Choice of a Crop Rotation 
Plan," Activity Analysis of Production and Allocation, Cowles 
Commission Monograph 13; New York, John Wiley & Sons, 1951, pp. 177- 
188. 

5. Dantzig, G.B., "Application of the Simplex Method to a Transporta- 
tion Problem," (ibid.), pp. 359-373. 

6, Melden, M.E., "Operetions Reseerch," Factory Management & Mainten- 
ance, October 1953, vol. 111, No. 10, pp. 113-120. 





Week’, In addition, the success of Operations Research and the atten- 7 
Ee à 

cent publicity has helped to focus attention to linear programming 

because this mathemetical technique is a useful tool for Operations 


Research teams. 








7. "Top Management by Methemetics," Business week, May 30, 1953, 
No. 1239, pp. 86-92, 





PURPOSE 

The majority of the literature to date has been slanted to the 
field of economics and mathematics and hence requires a certain degree 
of mathematical maturity for its understanding. A few articles 
eppearing in periodicals have intimated the possibilities of applica- 
tion of linear programming to Industrial Engineering problems. Between 
the two approaches lies a wide gap. 

The purpose of this thesis then is to attempt to bridge this gap 
by explaining what linear programming is, how the linear programming 
problem is interpreted mathematically, how to solve a simple problen, 
what the assumptions, limitations, and possibilities, of this technioue 
are, and how to set up for solution a few examples of feasible problems. 
The present treatment will endesvor to cover the material without 


requiring any mathematics beyond college algebra. 
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PRELIMINARY DISCUSSION 

Linear programming is a mathematical technique which may be 
utilized in the optimal planning of activities which are interdepen- 
dent. Programming in the title indicates the area of interest, the 
preparation of plans, which, when executed, will yield the best 
results. Programming is apparent in such areas, for instance, as 
mixtures of chemicals, loading of machine tools, assignment of person- 
nel, and shipping allocation. 

By linear is meant the restriction that the variables or 
unknowns must occur to the first power. No squares, cubes, etc., are 
permissable, nor may one variable be multiplied by another. Best or 
optimal results may refer, for example, to maximizing profits, or 
production, or to minimizing cost of production or storage. 

The linear programming technique is carried out as follows: 

1. A mathemetical model (or set of equations) is 
formulated from the word problem, this set of 
equations to follow a certain forn, 

2. The mathematical model is solved,usually by the 
simplex method. 

In regard to the formulation of the mathematical model it need 
not apply that the equations in the model express rigorously every 
facet, every fine point which could conceivably affect the problem 
at hand. Assumptions end approximations may be necessary. "Standard 
time", a statistical approximation, for example, may be employed in 
linear programming. Simplification of the problem may be advantageous, 
if not necessary. A certain phase of a production activity, for exam- 


ple, may be ignored if its effect on the total is small. When fornu- 





lating the model, therefore, assumptions, approximations and simpli- 
fications may be made provided the model retains enough predictive 
value to be useful. 

In the construction of the mathematical model there are two 
methods of approach. In the majority of the literature which, as 


has been stated, has evolved from an economics background where the 





interest lay in covering 8n entire economy, large formalized models 
of economic activities have ολα σα... This actívity analysis ap- 
proach, as it is called, has been extended to the theory of the firm 
so that many concrete and every day problems of industry may be 
investigated?. 

The second method of approach to the construction of the model, 
which is the one advanced by Charnes, Cooper and Henderson”, is "to 
approach each problem squarely, forming the model and calculational 
procedures naturally from the demands of the material." They feel that 
the activity analysis approach is awkward and artificial and "carries 
with it the danger of forcing the problem or overlooking important 
aspects which do not fit easily into prepared rubrics." 

The construction of a mathematical model, especially for a 
conplex problen, however, is an art and requires a keen mind coupled 
with a good mathematical and technological background. Examples of 


models will be presented later. 


8. Leontief, W.W., The Structure of the American Economy 1919-1939, 
second ed.; New York, Oxford University Press, 1951. 

9. Dorfmen, R., Application of Linear Programming to the Theory of the 
Firm; Berkeley; University of California Press, 1951. 

10. Charnes, A., Cooper, W.W., and Henderson, A., An Introduction to 


Linear Programming; New York, John Wiley & Sons, 1953, Preface. 








MATHEMATICS OF LINEAR PROGRAMMING 
Before illustrating the application of linear programming to problems 
it is necessary that the reader understand certain basic material of 
vector geometry and matrix algebra. Only the minimum of information 
necessary to the solution of linear programming problems will be 
presented. Mathematical proofs and amplifying material can be found 
in the material listed in the bibliography, especially in Part II of 
An Introduction to Linear Programming by Charnes, Cooper and Hender- 
gon, H 
Points and Vectors 
The representation of a point in a plane by (x, , y1) where 


X} is the x coordinate (distance from the origin or zero point along 


the x axis), and y, is the coordinate in the y direction, is well 
known. Just as (x), Was (xo, Yo) locate points in a plane, 
(x, , b 2.) and (χο, Yo, 2) locate points in the 3-dimensional space. 
(χι, Ji» 21)» in addition to representing a point, may represent a 
vector which can be visualized as an arrow whose tail is at the origin 
and whose head is at the point (x, Yi» 2i). 

Similarly, though less easily visualized, (x, ; T X35 re Xg) 
represents a point or vector in an n-dimensional space, where X4, X», 


through x, are real numbers called the coordinates in the dimensions 


1, 2, through n. Point or vector x may be represented either as 


(x, X555 X4) which is called a row vector, 


= 
aD 


or as , Which is called a column vector. 


EL, Ibid. 








The subscript under x designates the coordinate meant. A 
suverscript above identifies to which point one refers. Thus x and 
x^ means point 1 and point 2 respectively while x, means the second 
coordinate of point 3. 

Two or more vectors in the same n-dimensional space may be 


added by adding corresponding coordinates to form a new vector. Thus, 


η 1 & 2 Zu 2 
M x z ος x5, x3) and x = (x], Νο, x5) then x^, the sum of x* and 


2 
x is presented as X^» oc x, x5) = (1 4 x. x; + zc x, 4 x5) with 


1 : 


l ER, 1 Du l 2 
A AX 
This addition of vectors may be illustrated in ? dimensions as 
follows: Given: xl = (2,1); re (1, 3). Find x^, the sum of xl and 


x^. 


x! xt + e 


2:0]: 2 
{η . x x, + x^) 


το 1 + 3) 
= (3584). 


The point x63) may be found also by laying off from point x (2) a 


(2) 


vector parallel to and the same length as vector x' ', as shown in 


Figure ]. 








= 


4 x? = xt + x 
(3,4) 
2 
x 
2 (1,3) 
2 


(2,1) 


Fig. 1 Addition of Two Vectors 


À point or vector may be multiplied by a real number. Thus if 
x = (2, 1), then x= [(2)(2) , (2)(2)] = (4, 3). Such a miti- 
plication of a vector by a real number gives a new vector whose direc- 
tion is exactly the same as the original vector but whose length is 
equal to the length of the original vector multiplied by the nunber. 

A unit vector is a vector of value unity (one) in one axis only, 


and having value zero in every other axis. They are designated by 


ΠΡ 
9 , el ), d gn) etc., where the superscript designates which 


coordinate has value one. Thus in 2-dimensional space: οί) = (1l, 0) 


(2) 


and e = (0, 1) as shown in Figure 2. 











0 1 X 
Fig. 2 Unit Vectors in 2-Dimensional Space 
Similarly git ín an n-dimensional space is represented as: 
B. (0, 0,.., 0, 1, 0,.., O) where the l occurs in the ith dimension, 
Note that every vector in a 2-dimensional space (the plane) can 
be formed by a linear combination of the two unit vectors e) and 
e(2), For example: the vector x (1) = (3, 2) can be obtained from 
04 (003 (1,0) 42 (0,1) = (3, 2. 
Similarly, every vector in an s-dimensional space can be formed 
by a suitable linear combination of the s unit vectors. 
These s unit vectors in an s-dimensional space together form 


what is known as a linearly independent set. This means that for a 


linear combination of the s unit vectors, ay e (2) ta, a2) νο κ a, oU 
to equal the zero vector (0, C,........ , 0), the coefficients of the 


e's must each eoual zero. In 2 dimensions for example: 


84 ol) + a, eU) = "On. ©) 


ay (1, 0) + Bo (0, ντ, 9) 


OE (O ma) 20) ο) 





ο, τρια, p. L8, 
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= (a), ao) = (0, 0) only if a, - O and ao = 0. οἴ) and o) 
are thus a linearly independent set. 
No vector ín a linearly independent set can be formed from a linear 
combination of the other vectors in the set. For example, no linear 


combination of e) and οἵ2) can yield eU). 


This set of unit vectors is not the only linearly independent 
set. For example, in 2-dimensional space, (23, 1) and (2, 1) are 
linearly indepencent, for if a(3, 1) 4 b(2, 1) w (0, 0), then 

l. 3a + 2b = O and 

2.a+b=0 

From (2), b = -a, and substituting in (1), 3a - 2a =0 or 
a = O and thus b =0, 

Note, however, that in 2-dimensional space, (1, 3) and (2, 6) 
are not linearly independent but are linearly dependent. Each vector 
is a real number multiple of the other: 2(1, 3) = (2, 6). 

An additional consequence of linearly independent sets is that 


if a point P = a xl & a do... ...t a x? ( where xt, κά x is 


2 
a linearly independent set), then the coefficients 83, Boyer » 8. 
are the only coefficients of the xo 19 in terms of which one can express 


i 


P as a sum of multiples of the x Is, i.e., as a linear combination, 


The maximum number of points (vectors) in a linearly independent 
set in 2-dimensional space is 2, in 3-dimensional space is 3, and in 
m-dimensional space is m. In an m-dimensional space, any linearly 
independent set of m points is called a basis. There can be any number 


of bases but each basis of an m-dimensional space must consist of 





exactly m linearly independent points, 


Linear Transformations 

An important concept in linear programming is that of linear 
transformations. A transformation is a mapping or transfer of points 
from an n-dimensional space into a space of possibly different dimen- 
sions, each point assuming its new position in accordance with the 
directions specified by the particular transformation. For example, 
if the point x = (1, 3) were transformed by T to the point (2, 6) the 
transformation T of x ecuals (2) x indicated by T (x) s 2x. 

A linear transformation is merely a restricted tyoe of trans- 
formation and must conform to the following two requirements: 13 

l. T(ax) must equal aT(x). For example, if x = (1, 3), a = 2, 

and T specifies a multiplication by 3 


Tlax) = (3)(2)(1, 3) = (6, 18) 


aT(x) = (2)(3)(1, 3) = (6, 18) T thus conforms to 
the first requirement. 
2. T(x * w) must ecual T(x) + T(w). 
een. Specifies a multiplication by 3, 
x 2 (1, 3) and w = (2, 4), then 
T(x tw) = 3(1 +2, 3 #4) = (9, 21) 
T(x) = 3(1, 3) = (3, 9) and 


T(w) = 3(2, 4) = (6, 12) adding, 


T(x) $T(w) = (3 + 6, 9 #12) = (9, 21) 


T thus conforms to the second requirement. 





09. 1610, poA, 








Because the set of m unit vectors in an m-dimensional space is 
a linearly independent set, or basis, every point in the m-dimensional 
space can be expressed as a linear combination of the m unit vectors. 
Thus, if x = (3, 0, ?), x can be expressed as 


DS cM a le 


Χ 


3 
X= Cy et) where C, B Cy 5 ο C3 = 2 
isl 


It follows therefore that if it is known what the linear trans- 
formation does to the unit vectors, the linear transformation of any 
point is known. 

TE e) = (0,..., 1, 0,.., O) where the 1 occurs in the ¡th 
dimension, and i 2 l, 2,....., n, are unit vectors in an n-dimensional 
edie |) = (0,..,1, 0,....5°0), with 2 1, 2,:.., A, are UID 


vectors in the m-dimensional space M, then the linear transformation T 


m 
of οἱ 1) takes the form T οἱ 1) Ξ ai; gi. Thus: 14 
=] 


r e τα εἰ) κα el +0. ter y 
Te = i git) Bor g^ vue. A gm) zy 
ος) τα 6 tase) +. ba εἰ) ς γίν) 
T em) = a] g(t) + 82 £V) T NS añ gu) = y (a) 


The coefficients of the ES above when written by themselves 


can be wrítten in an array as follows: 


M e Ἱσ]αο- Po 4l. 
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457 822 ο . p 
Tac 

811 8.2 : A a 

a 

zd ο . e ad 


Such an array is customarily called a matrix of n rows and m colunns. 
The matrix is said to be associated (—~ ) with the linear transform- 
ation T. The entries in the us row are the coordinates of the point 
into which e (X) is carried by the transformation. Thus, e (2) which is 
the point (1, ©, 0,..., O) in the n-dimensional space N is transformed 
into the point y? whose coordinates are (241, 8125 **) 87? in the 


m-dimensional space M. 


Example: 
N is 8 3-dimensional space with 
e = (0, 1,0): 6 2 (0, 0088 
M is a 2<imensional space with 
ga) = (1, 0); τς = (0, 1), and: 
2 it 
Dr x 2 
2 3 
Then T οἴ) = 2 g(t) + 1 τ - y. a» (2,1) 
T a =] Sn 4352 g(?) k; = (0,2 
T εἰ”) -2 ga) t9 gl) = Tw» (2, 3) 


Figure 2 illustrates these transformations from space N into 


space M. 
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Fig. 3 Geometric Representation of Linear Transformations 


Since it is known what the linear transformation T does to the 
unit vectors in N, what T does to each point x is also known. If 
25 (χη, X2» X4) * 0,0, 1), 


D cape Tay x. Tet ναι το) 


moe ο ο CE. OL (2,09) 


G,2* 6» *[0«2,a«2] 20,0 


Points x and T (x) are illustrated in Figure 3. 


Convex Sets 


A convex set is a collection of points such that, if x and y 


are any two points in the collection, every point on the line segment 
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joining x and y is also ín the collection. Any point in a convex set 
‘which does not lie on the line segment joining some two other points f 
in the set is known as an extreme point. (See Figure 4). A rectangle, 
for instance, ís a convex set and the four corners are extreme points. 


A convex set which has a finite number of extreme points is given the 


special name of convex polyhedron. 


Extreme 
Point 


Convex Set Non-Convex Set 


Fig. 4 Convex and Non-Convex Sets 


In Figure 5 the convex polyhedron P4 Ρο P3 P is said to gen- 


erate the cone 0-P, P, P4 Pie This means that those rays emanating 
from the origin O which pass through the perimeter of the convex 


polyhedron 51 P form the sides of a cone. Such a cone is re- 


2 P3 P, 


ferred to as a convex polyhedral cone. 
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Fig. 5 Convex Polyhedral Cone 


Clearly the unit vectors in an m-dinensional space generate a 
convex polyhedral cone which will encompass all positive vectors in 
the m-dimensional space. 

The linear transformation of points has been discussed. Just 
as points in an n-dimensional space cen be transformed into image 
points in an m-dimensional space, so can a convex polyhedron. Indeed, 
the image under (resulting from) a linear transformation of a convex 
polyhedron is a convex polyhedron. An illustration of this is shown 
in Figure 2. The convex polyhedron in space N formed by the points 
ed, ea and eO) (a triangle) has been transformed into a triangle 


1l 
formed by the points y ^ y. and ya) in space M. 


Simplex Method of Solution 
The preceding mathematical discussion has presented some basic 
information necessary to the understanding of what a linear programming 


problem is and how it. can be solved. A simple numerical example will 
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now be presented and its solution will be obtained by the simplex 
method. 

The example consists of two inequations in three unknowns and 
a linear function (sometimes called a linear functional) of the three 
unknowns. The linear functional is to be maximized. A restriction is 
imposed that the unknowns when found must be non-negative, that is, 


must equal or be greater than zero (20). 


αι ..: οσο 
x *2y * 32 414 
f (x, y, 2) = x 4 y + 2z (maximum) 
wie x =O; y20: 220 
To convert the inequalities to equations a non-negative variable 
is added to each inequation and the model appears as follows: 


u LEM ca ct y 2 W710 
ie 


v + x +2y + ον. ον 


yv 


E 5, N 0 
3. f(x, y, 2) = x+y*+z + (0) u + (0) v to be a maximum. 

Stated in words the problem is to find the point or points in a 
5-dimensional space U which have non-negative coordinates (u, V, X, y,Z) 
satisfying equations (1) above and such that the linear functional 
(x + y +2) is a maximum. It can be shown that the point or set of 
points so described forms a convex set. The coefficients associated 
with the variables in equations (1) plus the right side of the equa-~ 


tions yield 6 column vectors P, in the 2-dimensional space V, Thus, 


i 
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l ‘0 12 [1 | 
er. E _ 12 10 
Pi z o | ; P, - a > P3 -= l| ; P, dd Ες > 3, ; and Pp = 14 
Equations (1) can be restated: 
Pz P, P3 Pa E Po 
x! ο) 2 1 2 10! 
E O A ro, πλ, 
or ur, 4 VP. + xP; + pe t 2P, - Po 


Recall the discussion of linear transformations. The coefficients 


in the above ecuations form a matrix associated with the linear trans- 





formation which takes the point Q = (u, v, x, y, z) in 5-dimensional 
space, U, into the point Ph = (10, 14) in 2-dimensional space, με” 
In the equations above, u is the coordinate along axis l, v 


along axis 2, in the space U, 


Poe) PA (1) g) 4 (0) gf?) 
Toe e (0) εἴ) 4 (2) gU) 
T e aep, αι (3) εἴ) (1) dz? 
TO - P = (1) g y aye 
T οί”) = Pe = (2) z Ὁ) + (3) zu where the 


coefficients of the g's are the elements in the re- 


spective vectors Fi 


There is more than one point Q in U having non-negative co- 
ordinates u, v, x, y, z, which will satisfy the equations, i.e., 
satisfy T (Q) = Po» and it can be shown thet these values of Y form 
Q = (u, vf x, y, Ὁ) = 


& convex set in U. For example: 


I Ibid. p.955. 
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ETC 1/3): a0, l, 0, 0, 0); and Q = (0, C, 3/2, 0, 9/2), 
are each answers which will satisfy the two eauations, as substitution 
back into the equations will show. All possible Q's which have non- 
negative coordinates u, v, X, y, and 2 form a convex set which, because 
it has a finite number of extreme points is a convex polyhedron. 

An important theorem states that a linear functional f (u, v, x, 
y, 2) defined on a convex polyhedron takes on its maximum (or minimun) 
at an extreme point of the convex set, 16 If the extreme points of the 
convex polyhedron generated by all values of Q which have non-negative 
coordinates could be found, the maximum value of the functional could 
also be found. That is to say, having the coordinetes of an extreme 
point, and substituting them into the functional, a value of the func- 
tional would be obtained. Comparing this value of the functional with 
the value obtained from every other extreme point of the convex set, 
one could pick out the extreme point whose coordinates yielded the max- 
imum value of the functional. These coordinates, or values for x, y, 
and 2 (u end v are of no interest, having been added merely to chanre 
the inecuations to equations) are then the reauired values for the 
optímum solution of the problem, 

The procedure, then, is to locate an extreme point, calcule te 
the value of its functional, proceed to another extreme point, anc so 
on, until the reouired maximum value of the functional is obtained. 

The extreme point which gives the maximum value of the functional is 


the solution. 


ντος, PF. 5». 








Recalling the discussion of linear independence, 1t was stated 
that in 2-dimensional space the maximum number of points or vectors 
in a linearly independent set, or basis, is two. The above points P 
are all in 2=cimensional space and subjecting any pair of points, P, to 
the criterion for linearly independent points as described on page $. 
it can be seen that any combination of two such points P is a linearly 
independent set. Any point, therefore, can be expressed es a linear 
combination of any two other such points P. For example, if x, Y» and 
z were each equal to zero, the ecustion could be stated as 

UI) “ο αυ 


HOR + ORTA OFP 


3 ? 5 = P 


2 0 


oruP, + vP, = Pg. Inserting the coordinates of 


2? and Po 


P, P 
u (1) + v (0) = 10 or ue ld 
u(0) + v(l) = Lor ve], 


The point Q m (10, y x» py 


3 ? 


0) is an extreme point. This 
point Q fulfills the requirements of the eauations of restrictions, but 
the value of the functional resulting from this point Q is not a max- 
— but is, instead, zero. The simplex method is designed to permit 
movement from one extreme point to another and in such a direction that 
an increase in the functional will result. 

To demonstrate the simplex method for the solution of this 


7 


ts are placed in a box arrangement. The Pa has been 


problem the P 0 


i 
moved into the first position. 








EI 


Table 1 





The original equations of the problem can be derived from the 
. box arrangement as follows: The single vertical lines represent plus 
signs and the double lines separating Po from the rest of the arrange- 
ment represent equal signs. Visualizing the elements of each horizon- 
tal row as being multiplied by their respective unknowns and including 
the equal and plus signs one can recreste the two equations. Thus from 
now 1: 

10 = u(1) + v(0) + x(2) + y(1) 4 2(2) 
and from row ?: 

w= u(0) #3 w(1 + x(1) 4 y(2) 1 2(3) 

An initial solution was obtained utilizing the basis set Py, Po. 
To show this, Table 1 is expanded into Table 2, where P, and Po have 


been placed to the left of Po in the stub column. 





Table 2 has the property of expressing any of the vectors in 


terms of the set Py and P5. Any vector as designated by the column 
headings can be stated by multiplying each element in the column by 








the vector on its row in the stub and adcins. Thus for Po the element 
on the first line under Po» 10, is multiplied by P] and added to the 
product of ihe second element, l4, with F>, 


Thus Po = (10) Py + (14) P, 
Similarly P) æ (1) EC (0) P5 
Po = (0) Py + (1) P, 
Pa = (2) P} + (1) Po 
A (A 
P (2) Py + (3) P3 


This follows from the fact that any point in a 2-dimensional 
space can be expressed as a linear combination of the points ina 
linearly independent set. Any of the above can be verified by inserting 


the numerical values of the P,'s. Thus for n 


4 = (2) U 4 alt) 


The numbers in the top row and far left column designated by c 


are the coefficients of the corresponding variables in the functional 


whose maximum is sought. Thus the number, 1, above P , and P. 16 
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the coefficient in the functional of x, y, and 2, resvectively. The 

coefficient of u and v is O as shown by the O above and to the left of 
Table 2 completes the first set of calculations. All of the six 

vectors (Pi, P», Pao Pao Pes and Po) appearing at the top of the table 


have been solved or stated in terms of the basis Py > Ρο and a feasible 





solution is indicated by the elements under Pos jon u e910, vie 14. 
Because Pas P, and Ps do not appear in the stub their corresponding 
variables x, y, and 2 are each equal to zero. The value of the func- 
tional at this stage is obtained by multiplying the feasible solution 
listed under Po by the functional coefficients listed to the left under 
c. The value of the functional is thus (10)(0) + (14)(0) = ©. 

Às was Stated, it is desired to move to another extreme poínt to 
investigate the value of the functional... The simplex method provides 
an algorism which yields a new value of the functional that is no smal- _ 
ler than the previous value. That is to say, the new value of the < | 
functional will either be larger than the previous value or possibly 
the same value as the previous value of the functional but will never 
be smaller. The new extreme point found will increase the value of one 
of the variables x, y, or z at the expense of cecreasing the present 
values of u and v namely 10 and l4, respectively. To find which vari- 
able to increase (what vector P3, Σι or Ρε to bring into the basis), 
the unit functional weights of x, y, and z are computed and compared. 
The largest value will designate the vector to bring into the basis. 

For exsmple, an increase in the size of x by one unit (from O to 
1) will increase the size of the functional by 1, the functional coef- 
ficient of x. But from Table 2, 1P3 œ 2 P} + 1l P3 oran increase 
in x by 1 will decrease u by 2 and v by 1. The decrease in the func- 
tional caused by reduction of u and v is equal to (0)(2) + (0)(1) = 0 
because the functional coefficients of u and v are each zero. A unit 


increase of x, therefore, will cause the functional to increase in 


value by 1. 
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Similarly, a unit increase of y will cause the functional value 


to increase by 


ea O υπ 
The remaining variable z, if increased by l, will likewise 
increase the functional by 1. 
e, È Cu u Cy V 
(1) (1) - (0) (2) - (0) (3) = 1 
Since the unit functional weights of x, y, and z are equal, the 
value of the functional will be increased the same amount no matter 
which variable is incressed (by inclusion of its P in the stub basis), 
hence it is immaterial which one is chosen. If the unit functional 
weights are not ecual, the one with the largest positive value desig- 
nates the F to include in the basis. 
Table < provides a means for a systematic computation of the 
unit functional weights as will now be shown. The table is reproduced 


here for convenience as Table 3. 





The values in the row designated by F are obtained as follows: 


each element in a column is multiplied by the c which appears to the 
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far left in the element's row and the products are added to give the 
value of F for that column. Thus the value of F under Po is equal to 
(10)(0) + (14)(0) = 0; F under Pz is equal to (2)(0) 4 (1)(0) «ο, 
etc. The F value under Py is the value of the functional obtained from 
the solution Q = (10, 12, δ, δ, 0). The F value in each of the other 
columns is the value of the decrease of the functional which would re- 
sult from the decrease in value of the two veriables u and v if a unit 
increase in the variable corresponding to that column were made. Thus, 
if x were increased from O to 1, the value of F, = 0 indicates that the 
reductions of u and v would have no effect upon the value of the func- 
tional. 

The elements in the row C-F are the unit functional woights, and 
are obtained by subtracting the F value in each column from the corres- 
ponding C value at the top of the column. 

The unit functional weights (C-F) are examined. The largest 
positive value determines which vector P to introduce into the busis in 
the stub to replace one P which is now in the stub. This replacement 
of an existing P in the basis by a new one corresponds to moving from 
one extreme point to another. By introducing the P which has the 
largest unit functional weicht (C-F) the largest increase in the value 
of the functional is insured. As long as there are positive values of 
(C-F) an increase in the functional is possible. 

In the present problem there are three ecual positive values of 
(C-F), each equal to 1. It is immaterial which one is chosen though 
the length of the computations might be shortened somewhat by a judi- 


cious choice. For illustration, Ps is chosen to replace one of the 
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vectors, Pz or Pos in the stub. 
The choice of which vector Pi ΟΥ̓ Ρο in tke stub to be replaced 
by the incoming vector Ps is made as follows: Divide each positive 


element in the Fs column into the corresponding element in the Po colunn 


and compare the results. Thus 


I5 
NIO 


4 1/3 


s 


The smallest result, here 4 1/2, designetes the vector Po of the 
stub as the one to be replaced by Pe. 4 1/3 or (14/3) is the new value 
of 2. Since the vector Ρο is being replacec, its varieble v decreases 
in value from 14 to 0. A blank table similar to Table 3 is now laid 


out but with P» in tke stub replaced by Ps. Note that the c æ l of Ps 


is placed to the left of Ps in the stub of Table 4. 


Table 4 





To complete the table, one recalls that a property of this ar- 
rangerent is that any vector as listed at the head of a column can be 
expressed as a linear combination of the two vectors in the stub. If 
the vectors P), Po, P3, Py, Pe, can be expressed as linear combinations 
of P4 and Pos the coefficients will be the elements of the table. From 


Table 3 the following equations hold: 





(a) Po = O D Pr 
(eer =e tere O PS 
ic) mo ορ T 1f 
5; rir 
.,, τ. 5», 
τ; 5. 
From equation (f) Po = -2/3 δι 4 1/3 Ps thus -2/3 and 1/3 are 


the elements in the new table under P.. Substituting -2/3 P3 * 1/3 Ες 


into equation (2) 


Po = 10P, 4 14 (-2/3 ss 1/3 Pe) 
or Po = 2/3 ΤΊ τ = Ες. Thus 2/3 and 14/3 are the elements 
to be filled in under Pos 
Similarly, 

Pi - l Pi + 0 Ες 

P3 4/3} 1 1/3 Ps 

P, = -1/3 P, + 2/3 Ps 

Po = 0 Pz + 1 Ps determine the elements for 


the other columns. The values of F and C-F are determined in the same 
manner as employed in Table 3. The completed table, Table 5 then 


appears as follows: 





Table 5 





Table 5 shows (in the Pg column) that a new solution u = 2/3, 
z = 14/3, with x = 0, y *0, v = 0, has been obtained. The F value 
below Po shows that the value of the functional has been increased 
from O to 14/3, and the two positive values in the C-F row show that 
further improvement is possible. The 2/3 value of C-F in column P4 
shows that P4 is the most profitable vector to introduce into the 
besis. As before the division of the elements under Po by those 
under Pz shows that Py is the vector in the besis to be replaced by 
P4. Carrying out the techniques employec previously Table 6 is 


constructed using the equations obtained from Table 5 to compute the 


elements in Table 6. 





Table 6 





Examination of Table 6 reveals a new feasible solution for Q, 
namely, Q * (0, 0, 1/2, 0, 9/2), and an increased value of the function- 


al, namely 5. The positive value of C-F, the 1/2 in column P,, indi- 


4 
cates thet a larger value of the functional can be obtained by introdu- 
cing P, into the basis. Because 3/4 in row P; is the only positive 
value in the P, column, Ρε is. the vector in the basis to be replaced. 
Carrying through the same procedure as before by computing the elements 
for the next table from the ecuations incicated in Table 6, the final 
table, Table 7 is completed 


Table 7 


y Z 
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The non-appearance of any positive values of C-F, the unit 
functional weights, indicates that no further increase in the value 
of the functional is possible. Reading in the Po column it is seen 
that the optimum solution is x = 2, y = 6, u #0, v "0,2 ΞΟ and the 
maximum value possible of the functional is 8. 

The procedure for arriving at the table elements, nemely by 
solution of equations from the prececing table, becomes extremely labor- 
ious when the number of variables and the number of equations 15 in- 
creased. A method exists that permits easy calculation of the elements 
by a systematic mechenical method. This method permits the solution of 
extremely large problems involving many variables and many equations in 
the original problem, 

The sample problem which has been sclved originated as a set of 
two inecuations. By introducing the additional unknowns to convert the 
inequations to equations the vectors Py «(1 4nd Ρο =(5) were intro- 
duced. These two slack vectors, as they are called, provided a known 
basis from which the first feasible solution (u z 10, v = 14) was ob- 
tained. When the size of the problem is increased the introduction of 
such a simple basis beccmes of more importance. If the original set of 
restrictive conditions are equstions instead of lnequations, slack 
vectors, now known as artificial vectors, may still be added to provide 
a basis. The associated variarles, however, must be included in the 
functional to be maximized anc eack such variable must have a negative 
coefficient of sufficient size so as to insure that the variable will 
not appear in the final answer. 


The above sample problem was concerned with the maximization of 





A fcd , 


pt that successive extreme points are found by an alporism which 
ures that each new value of the functional thet is obtained will be 


arger than any previous one. 








ASSUMPTIONS, POSSIBILITIES AND LIMITATIONS 

The bibliography contains accounts of the application of linear 
programming to industrial problems. In the main, however, the practi- 
cality of the technique suffers from certain restrictive assumptions, 
difficulties of model building, lack of proper data, necessity for 
complex computing machines and the necessity that a linear programming 
program be set up, if not operatec, by highly trained personnel. The 
last two may impose severe financial restrictions. 

The assumptions that must be met in order to apply linear progrem- 
ming are as follows :17 


a 
l. The assumption of linearity. A linear process is one whose ' 


ratio of output to input is ποσο ον t^ Thus for an operating machine 
tool, if an input of one unit of raw material produces one unit of 
output in a unit of time, an input of 100 units must produce 100 units 
of output in 100 units of time. Such en assumption neglects such 
aspects as rejects, tool wear, delay in delivery of material, set-up 
time, industrial fatigue, or unauthorized absence of the operator, etc. 
such non-linearity of a production process is, of course, nothing new 
and confronts procuction control personnel no matter what method they 


employ. Because of the large number of computations made in linear 


prograrming, however, and the fact that the results or decisions may 


17. Dorfman, R., Application of Linear Programming to the Theory of the 
Firm; Berkeley, University of California Press, 1951, p. €€C. 

18, This ratio is a result of the technical aspects of the specific 
process or activity. In the machining of metal, for instance, 
specified speeds and feeds in addition to other specifications 
determine the process. If the speed or feed were changed a 
technically feasible process might result. Such a modification of 
@ process, however, is, in linear programming terms , & new and 
distinct process. 








affect a large segment of a firm's production program it becomes neces- 
sary to insure a high degree of accuracy. Morgenstern states: 
Linear programning, or any other similar utilization of great masses 
of economic data, cannot be expected to make decisive practical 
progress until there is satisfaction that the data warrant the im- 
plied extensive and costly numerical operations ‚19 
Other examples of non-linearity are the reduction of sales price 
per unit for volume sales and the rise in productivity resulting from 
minor technological improvements. 
<. The assumption of divisibility. Any process may be operated 
at any positive level. This assumes, for instance, that if one large 
press can stamp out 100 fenders in an hour then one half of a press can 
stamp out 50 fenders in an hour. 
3. The assumption of additivity: 
It is assumed that two or more processes can be used simultaneously, Aur 
within the limitations of available resources, and that if this is ^ 
done the quantities of the outputs and inputs will be the sums of the œ 
quantities which would result if the several processes were used 
individually. 
4. The assumption that the nunber of processes is finite. ^ Be- 
cause of the restricted field of interest of many firns the processes 
or activities are indeed finite. Some industries, however, such as 
agriculture, oil refining, or the chemical industries have an infinite 


range of choices. 


Melvin E. Salveson has listed some of the problems which can and 


19. Morgenstern, O., "The Accuracy of Economic Observations," Activity 
Analysis of rroduction and Allocation, Cowles Comm. for Research in 
Economics, Monograph 13; New York, John Wiley and Sons, 1951, p. 283. 

20. Dorfman, K., Application of Linear Programming to the Theory of the 
Firm; Berkeley, University of California Press, 1951, p. 81. 








cannot be solved by linear programming. 1 He states that this tech- 
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nique can be used: 


1. 


to determine the optimum shop load during any one time period or 
over several time periods. 

to determine the optimum amount of overtime to use on any item or 
machine tool. 

to determine the optimum mix of commodities to make in the shop or 
factory. 

to determine the amount of rerouting of work in the shop. 

to determine the optimum level of inventory (raw, in-process, and 
finished). 

to determine the optimum distribution of production of commodities 
between time periods. 


He states that linear programming cannot be used: 


to determine a schedule for a shop. 

to select the best seguence of production in order to minimize, 
say, set up time. 

to eliminate need for expeciters. 

to avoid all conflicts in production, such as temporary bottle- 
necks end hence, its progrem is not necessarily always achievable 
(unless it is constructed with overly liberal delay allowances). 
to give delivery dates more precisely than the length of its time 
periods, 


The reasons for the non=applicability of linear programming to 


the five items above are as follows: 


ties. 


(a) Item 1 and item 5. These refer to the scheduling of activi- 


Production or any other activity of a linear programming problem 


when considered in relation to time is consicered to begin at the start 


of a specified tíme period, continue during the time period, at a con- 


stant rate and to reach completion at the end of the time period. This 


activity time requirement is faced in machine or shop loading problems 


also but suitable assumptions or restríctions in the formulation of the 





21. Salveson, M.E., "Mathematical Models in Management Programming," 


Journal of Industrial Engineering, March 1954, vol. 5, No. 2, 
Dp. 20% 





mathematical model may allow an operationally meaningful solution. 
An example of such a restrictive assumption will be given later. 

(b) Item 2. The selection of the sequence of production is 
determinec by the technology of the production processes. Thus, for 
example, the manufacture of a certain item might require that a milling 
-operation precede a drilling operation. This sequence is part of the 
technological requirements for the manufacture of the item. In adci- 
tion , set up time is a non-linear activity or function of a process. 
Thus one hour of set up time may be required for an operation regard- 
less of the level of activity of the process. 

(c) Item 3 and item 4. These restrictions on the epplicability 
of linear programming to production result from the inability of any 
mathematical model to account fully for the emergencies of the every 
day world. If every exigency were known ahead of time it might pos- 
sibly be incorporated in the model. The possibility of forecasting 
the occurence of every conceivable perturbation or disturbance in the 
orderly activities of nature is, of course, extremely slim. Many 
conflicts in production, such as machine breakdown or non-delivery of 
raw material due to personnel error, may be compensated for by con- 
structing the model with very liberal delay allowances. Expediters 
are necessary to keep the actual production as close to the program as 


possible. 





MODELS 

Many mathematical models have been presented in linear program- 
ming literature and work in the area of model building is proceeding. 
To show some of the rance of possibilities a few models are presented 
here. 

‘achine Assignment” 

A number of operations are to be performed in each of a number 
of mutually exclusive operation types. A number of machines of various 
types are available. The productivity of each machine type in each of 
the various types of operations is know. 

Required: to assign the machines to operations such that the 
total productivity of all machines on all operations is a maximum. 


Let: 
m z the nunber of types of machines. 


aj = the number of machines of type (i) where 
lo eae, y ΜΟΝ 
n = the number of types of operations. 
b. = the number of operations of type (j) whore j = 
ee ee un. 
C4; = the productivity of a machine of type (i) in 
performing operation (j). 


the number of type (i) machines to be assigned to 


X14 
a type (3) operation, 


22. Harrison, J.O., Jr., Linear Programming and Operetions Research, 
Informal Seminar in Operations Research, Paper No. 2; Baltimore, 
Johns Rapkins Univ,, Octeber 13, 1953, pp. 17-13, 





The X14 are to be selected such that: 


2 
X14 ara) 


n 


X14. $ 8, 


3: 
lM 


) ES - 5 
DN n 
051 AS maximum 
1-1 J= 


Production Line Loading to Meet Sales Requirements” 

Given: one production line producing one type of product which 
sells for a fixed unit price. The unit costs of regular time produc- 
tion, overtime production, and storage are known. The rates of produc- 
tion per unit time are known. A sales demand during each of a number 
of successive time periods is known. The time period is one month. 

Desired: a production loading plan which will meet the sales 
demand and minimize the combined costs of production and storage. 
Inventory is taken at the end of each month. 

let: 

κ = number of time periods (months) to load. 

S; = number of units of finished product to be sold during the 

ith period. 
Ig = initial inventory. 
R = maximum mmber of units which can be produced during a 
month on regular time. 
U = maximum number of units which can be produced during a 
month on overtime. 


cost of storage of one unit of product for one month. 





23. Ibid. pp. l4-15. 





C, = unit production cost utilizing regular time, 


Cy = unit production cost utilizing overtime. 


Ping: , 
X; = number of units to be produced on regular time during 
the i?! month. 
ὅ = number of units to be produced on overtime during the 
11h month. 
(a) zn 9 Ji 2 0 (i = 15 21.00 als 
(b) X. € R 
(ο) σι £ U 
d a 
(a) cat ? Sa -Io 
Ez a=/ 
(e) (Cr Iy + Cy x, $ Cy y4) to be minimized, where 
1-71 i 
nt (Xa + Ya Sa 


a=/ 
A 


Equation (d) is obtained as follows: 

The inventory during any month, a, is denoted by I,. Since 
the inventory at the end of any month is equal to the inventory at the 
end of the previous month plus the monthly production minus the monthly 
sales, 

la = la] t Xa * Ya - Sa 

NOSE this equation from a = l to a * i yields equation (f). 
Since the right side of equation (f) must be equal to or greater than 
zero ( 20), or i i 

1) + (Xa 4 Ya) - Sa O 


=] =f 
A slight rearrangement yields ecuation (d). 





å Numerical Example - Product 
Assembly Over Consecutive Time Periods/^ 

To illustrate the complexity to which even problems of small 
size are subject a numerical example is presented. 

Assumptions: a factory making three products, each of which 
can be final products and two of which are in addition intermediate 
products. 

Product 1 (Fi) is an assembly of one unit each of P, and P3, 

Product 2 (P5) is composed of 2 units of Product 3 (Ea). 

There are two types of operations, type 1 (Τη), and type 2 (Το), 

One hour of T, is required for one unit of Pie 

One hour of T plus one hour of To is required for one unit of Po. 

Two hours of To are required for one unit of Ps. 


The above assumptions are summarized in the following table; 


To make To make To make 
one unit Pi one unit P, one unit P4 
P4 (unit) 
P» (unit) l 
Pa (unit) l 2 
Τη time (hours) 1 l 
T, time (hours) 1 2 





24. Jackson, J.P., 'Mathematical Models with Examples from Linear 
Progremning," Industrial Logistics Research Froject Report No. ὁ; 
Los Angeles, Univ, of California. 
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Further assumptions: 


3000 hours per week of T, are available. 


1 
¿000 hours per week of T, ere available. 
3000 units of P3 are in storage at the commencement of the 
production period. 


Net profit from the sale of: 


Py = $3 
Py = $5 
P oes 


The net increase in profit from the sale of: 
PR = l4 - 5 - 1 = $8 and of 


5 =- ο) = $3 


Itens produced during any given week cannot enter into further 
production until the next week. 16, units of Po made during the first 
week cannot be assembled into Pz until the second week, 


Recuirement: to operate for three weeks so as to maximize profits. 


Let: 
xt = production of Pz in week t: t = 1, 2, 3. 
y? = " " P,inweek t: t = 1, 2, 3 
+ 

Z’ m E "P3 inweek t; t = 1, 2, 3 

u* s number of units of P4 stored during week t; t = 1, 2, 3 
yt - " n " " Ρο " " "o" bes 1,2 408 
wt = " " n n P3 n " "oom, { = 1, 2,093 
at = number of hours of idle time of Ty in week t; t°@ 1, 2,9; 


b? = number of hours of idle time of Tə in week t; t = 1, 2, 3. 





From the table it can be seen that: 


O = number of units of P4 consumed in week t; Ὁ = 1, 2, 3. 
t 
x e " " " " Po " " p cc Ds 2 3, 
xt42yt = " " " " Pa " " "obo t E ], ΠΣ 





In addition: 


O = number of units of P, available at beginning of week 1. 


Bos nono nup : "oom "on oon 
EV o: oo ας ο o "oon ho om om, 

ut y xt = P. available at beginning of week (t * 1); t = 1, 2, 3. 

ve. yt = >, n " " " " (t $1; t - 1, EN 

πὶ 127 = PS " " " " " (t&1); te 1 2 


From the above restrictions the following equations can be obtained: 
= 0 
vi + yl = 0 


wl + xl 4 oyl = 3000 


ιτ. ut = ut t xt-l, per 3% 

yt + xt m yt-1 + yt t a = Ὧν 

wt t xt 4 2y* = imd t αὔτ] πα 2,3% 

Availability of machine time leads to the following equations: 
1  * + yt 4 at = 3000 ν Ἡ 1, δ, 3 


In addition: 
III. xv, yt, αὖ, αὖ, vt, wt, at, pius 0; towel E. 2,011 
Finally the net profit: 
IV. P= (xt + x7 4 x3) * 3(yl * y? 5 y?) + (zl & z? « z3) 
Thus, to operate for three weeks so as to maximize profits calls 


for maximizing IV subject to the restrictions imposed by I, II, and III. 


The preceding problem by J.R.Jackson can now be solved in the 





folloving way: 





Substituting the valves for (t) the following equations are 
obtained: 

ul = 0 

vi + yl z ο 

wl + xl + 21 = 3000 

u? = ul 4 xl 

uve = u? + x2 

v^ * x? - τὶ 4 yl 

v3 + 3 = v2 + y? 

We my ox* 4 2y^ = yl 4 zl 

wo & x3 4 2y3 m w^ 4 z^ 


xl + yl + al = 3000 


> 
N 
+ 
«4 
N 
+ 
f5 
i 


3000 

x + y? + ad = 3000 

yt + 222 + pl = 4000 

y2 + 222 + b* = ¿000 

y? $ 222 * b^ = 4000 

Since ul = O it may be eliminated, and since vi = -xl and u? = xt, 
vi and u“ may be replaced by -xl and xl respectively. The result is a 
set of 12 equations as follows: 

l. wl + xl + 23y} = 3000 

2.1} - χὶ - χ = 0 

3.2 ı κ2 -y amo 

4. v - ve = y% = 0 


2 we + x* 4 2y^ ER OT gt = 0 


6. w 4 3 à o - 2 - 22 = 0 


7, xl + yl + al 2 3000 
8. x + y2 & a2 = 3000 


3000 


9, x? 4 y2 4 a? 
10. yl 4 21 + bl = ¿000 
ll. y^ + 222 + b? = 4000 
12. y) + 2? & t) z 4000 

These 12 ecuetions are displayed in the femiliar arrangement of 
Table 8. It appears that it will be possible to rearrange the table 
so as to provide twelve linearly independent vectors to form the initial 
basis. The below listed steps will yield an equivalent matrix which 
will have the desired form. 

l. Replace row 6 by the sum of row 6 and row 5. 

2. Interchange columns v? and wl. 

3. Replace row 4 by the sum of row 4 and row 2. 

4. Replace row 5 by the sum of row 5 and row 1. 

5. Replace row 6 by the sum of row 6 and row 1. 

6. Interchenge columns v* and wl. 

7. Interchange columns u? and wl. 

Moving the basis to the left side of the table, adding the 
functional coefficients, the stub column, and the F and C - F rows 
completes the first set of celculations as shown in Table 9. Suc- 
ceecing tables can then be utilized to carry out the steps of the 


simplex method, 
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To show the amount of computation required to solve this problem 
by an automatic computer recall that advancing from one table to the 
next is in effect the movement from one extreme point to another. To 
go from one extreme point to another by the following steps requires 
approximately the number of computer operations as indicated for each 
step. 
l. the calculation of the unit functional weights (C - F). 
286 operations. 

2. the inspection anc comparison of the unit functional weights 
anc the selection of the one with the largest positive value. 
This determines which vector is to be placed in the basis 
36 operations. 

3. The selection of the basis vector which is to be repleced by 

the incoming vector. 60 operations. 

4. the calculation of the new value of the functional, 2 opera- 

tions 

5. The determination of the column vectors in terms of the new 

basis. 242 operations. 

The total number of operations required then to advance from one 
extreme point to the next is 626, 

Computational experience indicates that, in general, the number 
of extreme points which must be examinec is approximately twice the 
number of equations. In this problem then, the total number of computer 
operations is (2)(12)(626) z 15024. 

A Card Programmed Calculator (C.P,C.) which is a general purpose 


floating point digital computer that can perform 100 operations per 





minute would then take approximately two and one-half hours of actual 
computer operation to solve this problem, This computer, however, can- 
not store enough information to allow continuous operation. The manual 
bookkeeping required for the removal and re-insertion of information 
doubles or triples the time required, and this problem, therefore, 


would require between five and seven and one half hours for solution. 





48 


CONCLUSION 

Linear programming is thus a mathematical technique to find the 
values of non-negative variables or unknowns which are subject to linear 
restrictions and such that the values will optimize a linear functional. 
It has been shown that the linear programming problem can be visos lied 
as a geometric representation of vectors or points in imaginary spaces 
of many dimensions. The solution of a simple numerical example has 
been demonstrated using the simplex method. 

The requirements that must be met before this technique can be 
used for practical problems have been discussed. It can thus be seen 
that linear programming has e restricted usefulness for Industrial 


Engineering applications. As Joseph O. Harrison says 12 


The practical difficulties encountered in applying linear programming 
are threefold: (1) the expression of realistic objectives and con- 
straints in measurable terms, (2) the determination of suitable numer- 
ical values for coefficients, and (3) the computational labor required 
to numerically execute large linear programming problems. 

Like any mathemstical technique, linear programmáng must not be 
considered a panacea for all the problems of Industrial Engineering. It 
is restricted in the scope of its application both by the nature of the 
problem and by the economic considerations of computation. 

The history of the application of linear programming to practical 
problems covers too short a period of time and the application has been 
too limited to permit a complete evaluation of the usefulness of this 


technigue to Industrial Engineerinz. It can be expected, however, that 





25, Harrison, J.0., Jr., Linear Programming and Operations Research, 
Informal Seminar in Operations Research, Faper No. 2; Baltimore, 
Johns Hopkins Univ., October 13, 1953, p. 18. 
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the considerable amount of study now being devoted to both the theory 
and the practical aspects of linear programming will result in the 
overcoming of many of the difficulties aeo its usefulness, and 
that Industrial Engineering will make more and more use of this tech- 


nigue in the future. 
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